Accepted for publication in the Astrophysical Journal 
Preprint typeset using I4TgX style emulateapj v. 10/09/06 



THE NON-ISOTHERM ALITY AND EXTENT OF GALACTIC DIFFUSE HOT GAS TOWARD MRK 421 

Yangsen Yao 1 and Q. Daniel. Wang 2 

Accepted for publication in the Astrophysical Journal 



o 
o 

O 

Q 

oo 



ABSTRACT 

Diffuse hot gas can be traced effectively by its X-ray absorption and emission. We present a joint-analysis 
of these tracers to characterize the spatial and temperature distributions of the Galactic hot gas along the sight- 
line toward the nearby bright active galactic nucleus Mrk 421. We also complement this analysis with far-UV 
O vi absorption observations. We find that the observed absorption line strengths of O vn and O vm are 
inconsistent with the diffuse background emission line ratio of the same ions, if the gas is assumed to be 
isothermal in a collisional ionization equilibrium state. But all these lines as well as the diffuse 3/4-keV broad- 
band background intensity in the field can be fitted with a plasma with a power law temperature distribution. We 
show that this distribution can be derived from a hot gaseous disk model with the gas temperature and density 
decreasing exponentially with the vertical distance from the Galactic plane. The joint fit gives the exponential 
scale heights as ~ 1 .0 kpc and 1 .6 kpc and the middle plane values as 2.8 x 10 6 K and 2.4 x 10~ 3 cm" 3 for the 
temperature and density, respectively. These values are consistent with those inferred from X-ray observations 
of nearby edge-on galaxies similar to our own. 

Subject headings: Galaxy: disk — Galaxy: structure — X-rays: ISM — X-rays: individual (Mrk 421) 
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1. INTRODUCTION 

X-ray absorption lines of O vn, O vm, and/or Ne ix at 
z ~ have been detected unambiguously toward several bright 
AGNs, but the nature of these highly-ionized species remains 
uncertain. Because of the limited spectral resolution of the 
current grating instruments on board the Chandra and XMM- 
Newton X-ray observatories (FWHM ~ 50 mA, correspond- 
ing to a velocity dispersion of 750 km s" 1 at 20 A), however, 
the line width and centroid position do not directly constrain 
the absorbing gas to be any significantly better than a few 
Mpc. Some authors have attributed these z ~ absorptions 
to the so-called warm hot intergalactic medium (WHIM) in 
the Local Group or even beyond (e.g., Nicastro et al. 2002; 
Rasmussen et al. 2003; Williams et al. 2005). In this inter- 
pretation, the hot phase medium could contain large amounts 
of baryonic matter, accounting for the bulk of the "miss- 
ing baryons", as predicted in numerical simulations of large- 
scale structure formations (e.g., Cen & Ostriker 1999; Dave 
etal. 2001). 

But, the z~0O vii, O vm, and/or Ne ix absorption lines 
have also been observed in X-ray spectra of the LMC X-3 
at distance of 50 kpc (Wang et al. 2005) and many Galac- 
tic sources (e.g., Futamoto et al. 2004; Yao & Wang 2005; 
Juett et al. 2006). The strength of the absorption toward such 
a source, depending its Galactic coordinates and pathlength, 
can be substantially greater than those toward the AGNs. A 
crude characterization of the global distribution of the absorb- 
ing hot gas gives an effective scale length of several kpc, 
suggesting that much of the observ ed absorptions originates 
in the Galaxy (I Yao & Wang 1 120051) . The recent surveys of 
AGN spectra with the Chandra and XMM-Newton archive 
data show that t he O vn absorption line is present in all high 
quality spectra dMcKernan et al. 1120041: iFang et ani2006l) . in- 
dicating a covering fraction of O vn absorbing gas >63% in 
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the sky. This high covering fraction, together with an estima- 
tion of the total mass of the absorbing gas and with a compar- 
ison between emission and absorption measurements, further 
supports the scenario that t he X-ray absorbe rs are primarily 
associated with the Galaxy dFang et al. 1 12006). There is really 
little doubt as to the presence of the Galactic contribution to 
the X-ray absorption. The question is how much. 

While the line absorption depends on the ionic column den- 
sity and the dispersion velocity of the intervening hot gas, the 
emission is sensitive to the number density and temperature 
of the gas. Therefore a combination of absorption and emis- 
sion measurements naturally provides a cons traint to the scale 
length of t he absorbing/emitting plasma. iRasmussen et al. I 
(2003) and Futa moto et al. I (12004) have indeed compared the 
observed O vn absorption with the back ground O vn emission 
line intensity (iMcCammon et al.~1l2.Q02h to estimate the scale 
length of the hot gas around the G alaxy. Similar estimates are 
also made by IFang et al. I (120061) . But, so far such estimates 
all depend on various assumptions such as absorption in the 
linear regime and isothermality of the hot gas, which may not 
be valid (see 93. H and §4 for further discussion). 

Here we report the first systematic joint-analysis of X-ray 
absorption and emission data, complemented by a far-UV 
OVI absorption measurement, in the direction of Mrk 421. 
The joint-analysis allows for tests of the assumptions and pro- 
vide tight constraints on our spectral models. Our goals are to 
characterize both the scale length and global thermal property 
of the z ~ absorbing/emitting hot gas toward Mrk 421. 

Throughout the paper, the statistical errors are quoted 
at the 90% confidence level, unless being pointed out 
specifically. We al s o ado pt the solar abundance given by 
lAnders & Grevessel (119891) for oxygen 3 . 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra X-ray Absorption Data 



3 We note that Asplund et al. 1 2005) propo s ed a r evision of the commonly 
used solar abundances oflAnders & Grevesse I fT989T) . but this revision is still 
under debate I Antia & Basu 2006). We therefore still use the old values. 
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TABLE 1 

Chandra OBSERVATIONS 



ObsID 


Obs. Date 


Grating/Detector 


— 

Exposure 

(ks) 


Used 


457 


1999 Nov 5 


HETG-ACIS 


19.83 


No 


1714 


2000 May 29 


HETG-ACIS 


19.84 


No 


1715 


2000 May 29 


LETG-HRC 


10.89 


No 


4148 


2002 Oct 26 


LETG-ACIS 


96.84 


Yes 


4149 


2003 Jul 1 


LETG-HRC 


99.98 


Yes 


5171 


2004 Jul 3 


LETG-ACIS 


67.15 


Yes 


5318 


2004 May 6 


LETG-ACIS 


30.16 


Yes 


5331 


2004 Jul 12 


LETG-ACIS 


69.50 


Yes 


5332 


2004 Jul 14 


LETG-ACIS 


67.06 


Yes 



Mrk 421 is one of the brightest quasars at z = 0.03 (l,b = 
179.°83,65.°03), and its X-ray flux varies by up to ~ 50 times 
between quiescent to burst states. Table [T] lists the existing 
nine Chandra observations of Mrk 421, covering various flux 
levels of the source. Two observations used the High En- 
ergy Transmission Grating (HETG; Canizares et al. 2005) 
with the Advanced CCD Imaging Spectrometer (ACIS-S; 
Garmire et al. 2003) and the other seven used the Low En- 
ergy Transmission Grating (LETG) with either the ACIS-S 
(five times) or the High Resolution Camera (HRC-S; Mur- 
ray & Cappell 1986). The LETG provides - 30-40 cm 2 
effective area around O vn and O v m Ka waveleng t h with 
a m odest resolution of ~ 750 km s" 1 . Nicas tro et al. I (120051) 
and I Williams et all (12005b have reported three (ObsID 1715, 
4148, and 4149) of these observations. In this work, we com- 
bine six long exposures of these observations to enhance the 
counting statistics (by ~ 30%); the other three observations, 
which we do not use, have short exposures, contributing only 
< 5% to the total recorded counts around the wavelength of 
O vii Ka (the most significant absorption line). Furthermore, 
the LETG-ACIS combination suffers no grating order over- 
lapping complexity as the LETG-HRC does, allowing for a 
sensitive cross-checking of the measurements between the 
ACIS-S and HRC-S observations (see below an d §4). 

We follow the same procedure as described in lYao & Wang I 
(2006) to re-process the individual observations. For the 
ACIS-S observations, we extract the first grating order spectra 
and calculate the corresponding response matrix files (RMFs) 
and auxiliary response functions (ARFs). For the HRC-S ob- 
servation, we extract the total grating spectrum and calculate 
the RMFs and ARFs from the first to the sixth grating orders, 
and then use IDL routines from PINTof ALE package to com- 
bine all the RMFs and ARFs pairs to form an order-integrated 
response file accounting for the higher order overlapping in 
the spectral fitting. To improve the counting statistics, all the 
spectra, after a consistency check, are then co-added and the 
corresponding ARFs and RMFs are weight-averaged accord- 
ing to the exposure and global continuum intensity obtained 
from a model fitting to each spectrum. In the final co-added 
spectrum, we obtain about 2700 and 1600 counts per channel 
with a channel width of 12.5 mA around wavelength of 18.5, 
and 21.5 A, where our interested absorption lines are located. 

Our global continuum modeling is conducted in the 2- 
22.5 A range. We adopt a power law model with a fore- 
ground cool gas absorption to characterize th e overall spec- 
tral shape of Mkr 42 1 . A s in previous studies (Nicastr o et al. I 
120051: 1 Wang et al71 12005), in order to obtain an acceptable 

4 http://hea-www.harvard.edu/PINTolALE 



global fit, we include 11 broad (a > 1500 km s ) Gaussian 
profiles to account for spectral deviations from the adopted 
power law model and for various calibration uncertainties 
in detector gap/node regions. The Ne ix, O vm, O vn Ka, 
and O vn K/3 absorption lines are clearly visible in their 
corresponding rest frame wavelength (we do not use other 
marginally significant z ~ absorption lines that have been 
reported by Williams et al. 2005 and by Nicastro et al. 2005). 
These absorption lines are characterized with narrow Gaus- 
sian profiles in the continuum fit. We obtain a final spectral fit 
with x 2 /d.o. f = 1987/1577, which, given such high counting 
statistics, is reasonably acceptable and in particular is good 
enough for accounting for the higher order overlapping in 
the LETG-HRC data over the whole wavelength range. The 
X 2 /d.o.f in the 12.5-14.5 A 18.2-20.0 A, and 20.5-22.5 A 
ranges are 202/148, 180/148 and 202/147, respectively. 

2.2. Diffuse X-ray Emission Data 

We use the diffuse 3/4-keV band (R4+R5 bands) back- 
ground intensity from the ROSAT All Sky Survey (RASS; 
Snowden et al. 1997) to estimate the emission from the hot 
gas in the field of Mrk 421. The background at lower ener- 
gies (i.e., 1/4-keV band) is subject to strong absorption by the 
cool interstellar medium and therefore traces mostly the hot 
gas within the Local Bubble. The background in the 1.5-keV 
band is primarily point-like in origin. To avoid any signifi- 
cant contamination from the bright Mrk 421 due to the ex- 
tended point spread function wing of the ROSAT, we use an 
annulus between 0.°75 and l.°75 radii around the source to 
obtain a 3/4-keV background intensity estimate as ~ 1.06 x 
10~ 4 counts s" 1 armin" 2 . This background intensity contains 
a significant contribution from unresolved point-like sources. 
The extragalactic point-like source contribution can be char- 
acterized as a power law function with photon index 1 .4 and 
a normalization 10.9 photon s s" 1 keV" 1 cm" 2 sr" 1 at 1 keV 
dHickox & Markevitchl l2006). which is equivalent to 5.61 x 
10~ 5 counts s" 1 armin" 2 in the RASS 3/4-keV energy band 
(after considering the foreground cool gas absorption = 
1.43 x 10 20 cm" 2 ; Dickey & Lockman 1990). The Galactic 
component (due to stars) is ~ 6 x 10~ 6 counts s" 1 armin" 2 a t 
the Galactic latitude of Mrk 421 dKuntz & Snowden Il200ll) . 
Therefore, we estimate the background intensity due to the 
Galactic diffuse hot gas as 4.6 x 10~ 5 counts s" 1 armin" 2 . 
Here we have neglected a potential contribution from charge 
exchange interactions of solar wind ions with the local neutral 
interstellar medium (ISM; e.g., Pepino et al. 2004). 

The 3/4-keV background is predominantly due to the O vii 
and O vm Ka line emission. The best existing spectro- 
scopic measurement of such line emission is from the micro- 
calori meters flying on a sounding rocket dMcCammon et al."1 
120021) over a region of ~ 1 sr, centered at l,b = 90°, 61°, 
which is at nearly the same Galactic latitude as the Mkr 421 
direction and is certainly away from abnormal regions such 
as the Galactic Bulge, the North Polar Spur, and the Cygnus 
Loop. The measured line intensities are /ovii = 4.8 ± 0.8 and 
^oviii = 1.6 ±0.4 photons s" 1 cm" 2 sr" 1 (la errors). We as- 
sume these intensities are typical for the hot gas at high lati- 
tudes. The effects of the foreground cool gas absorption are 
expected to be small 5 and therefore have not been consid- 

5 The neutral hydrogen column density toward the rocket experimental 
direction is 1.9 X 10 20 cnT 2 . The absorption correction will increase the 
OVII and OVIII line intensities by 18% and 12% respectively; these changes 
are comparable to or less than the corresponding \a statistic uncertainties. 
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ered in this work. We use these two line intensities, imple- 
mented in a two data point spectrum, to characterize the hot 
gas emission. This spectrum will be jointly fitted with the 
absorption data to constrain the emitting/absorbing gas prop- 
erties, whereas the broadband RASS 3/4-keV intensity will 
be utilized as a consistency check. We also fake a response 
file containing two delta functions for this spectrum, so this 
emission-line-intensity spectrum can be directly compared to 
our emission line model in the joint analysis. 

2.3. FUSE O vi Absorption Data 

FUSE 6 observed Mrk 421 four times with 29 individual ex- 
posures and an actual integrated exposure of 84 ks in total, all 
configured with the target aperture LWRS (30" x 30"), pro- 
viding effective area of ~ 50 cm 2 with spectral resolution of 
~ 20 km s" 1 around the O vi doublet (1032 and 1038 A). 

We only use the data from the LiFla segment. We cali- 
brate the wavelength for each exposure with CalFuse v3.0 7 , 
and then repeatedly cross-correlate the wavelength scale of 
different exposures and co-add them to form a final com- 
bined spectrum using Don Lindler's IDL Tools 8 . In order to 
jointly-analyze the FUSE and X-ray spectra together in soft- 
ware package XSPEC, we transfer the FUSE FITS spectrum 
to a format that can be recognized in XSPEC. We adopt a 
Gaussian line spread function with 20 km s" 1 (FWHM) as the 
spectral response of the FUSE spectrum. 

We use a power law to approximate the continuum and a 
Gaussian profile to characterize the O vi absorption. We find 
that one broad Gaussian component is adequate to describe 
the total absorption. We obtain an EW of the O vi line that is 
consistent with the previously reported value (e.g., Wakker et 
al. 2003; Williams et al. 2005). As the line is well-resolved, 
we measure the dispersion velocity and O vi column den- 
sity of the absorbing gas as v_ 9 = 90(84, 97) km s" 1 and 
log[jV vi(cnr 2 )] = 14.44(14.41, 14.47). 

3. DATA ANALYSIS AND RESULTS 

Our analysis includes the absorption and emission lines of 
oxygen only to avoid uncertainties in the relative elemental 
abundances of the hot gas. We assume that the hot gas is in 
a collisional ionization equilibrium (CIE) state and that the 
X-ray absorption contains only one velocity component. In 
reality, each observed line could consist of multiple veloc- 
ity components. With the still limited spectral resolution of 
current X-ray data, we decide not to investigate this complex- 
ity and leave it to the future when better data become avail- 
able. The absorption lines of various transients are jointly 
fitted with our model absline 10 , implemented in the spectral 
analysis software XSPEC. 

3.1. Isothermal Model 

We first conduct a joint analysis of the O vn Ka, O vn K/3, 
and O vm Ka X-ray absorption lines (Fig.[T|). The goodness 
of the spectral fit with the absline model is comparable to that 
obtained with the Gaussian profiles Q2.1I ). From the fit, we 



6 For detail description of FUSE, please refer to Moos et al. 1 2000) 

7 http://fuse.pha.jhu.edu/analysis/calfuse.html 

8 http://fuse.pha.jhu.edu/analysis/fuse_idl_tools.html 

9 We use i_ instead of the more generally used b to denote the standard 
velocity dispersion to avoid the confusion with the Galactic latitude. 

10 For a detail description of the absline mod el and its application in mul - 
tipl e absorption line di agnostics, please refer to Yao & Wang 1 20051 l200rj| 
an d Wan g et all JlQOBT) . 



obtain the dispersion velocity v_ , O vn column density Novn, 
and the temperature T of the intervening gas. We can then 
infer the O vi, O vm, and hydrogen column densities (Novi, 
Novm, and Nh respectively) of the absorbing gas (Tabled. 



0VII K/3 18.629 A 




Velocity (km s ) 



FIG. 1 . — The joint fit of the X-ray absorption lines with an isothermal 
absline model. The bin size is 12.5 mA; the ARF and RMF have been applied 
to the model. 

We further jointly analyze the O vn emission line inten- 
sity (without including the O vm line; _2.2b with the X- 
ray absorption lines to estimate the scale length of the emit- 
ting/absorbing gas. The oxygen emission line intensity is 

/[ph cm -2 s _1 sr _1 ] = — / A(T)n e ndl, (1) 
4tt Jo 

where Aq is the oxygen abundance (in unit of the solar value 
and assumed to be 1 in the paper), A(T) is the line emissiv- 
ity , and L is the effective path length of the gas. Assuming 
that the hydrogen number density can be characterized as 



n = riQe 



(2) 



where z is the vertical distance away from the Galactic plane, 
«o and h n are the mid-plane value and the scale height, and 
£ is the volume filling factor that is assumed to be 1 in the 
paper, we can integrate Eq.Q]along a sight-line with a Galactic 
latitude b to get 

r 1.2A N 2 H A(T) 



8ttL 



(3) 



Here the factor 1.2 accounts for the helium contribution to the 
electron density, L = h n £/ sinb, Nh = «o^n_/ smb. We con- 
struct a simple emission line model based on Eq. [3] and use 
it to model the emission line intensity spectrum we have con- 
structed ( jj2.21 >. In the joint analysis, and T are the com- 
mon fitting parameters in both the absorption and emission 
line models, whereas L only appears in the emission model. 

11 The definition of A(r) is adopted f rom 

http://cxc.harvard.edu/atomdb/physics/physics_units/physics_units.html 



4 Yao & Wang 



TABLE 2 
Joint analysis results 



v b log(T) log(A'ovi) logCAfavn) log(A'ovm) log(Ar H A ) D 

(kms 4 ) (K) (cm 4 ) (cm 4 ) (cm 4 ) (cm 4 ) 

Isothermal Modeling 

abs! 64(48,104) 6.16(6.10,6.21) 13.51(13.32,13.70)* 16.00(15.82,16.17) 15.18(14.78,15.41) 19.14(18.98,19.30) 

Non-isothermal Modeling 

abs. + emi. 64(48, 103) 6.44(6.36, 6.54) c 14.76(14.53, 14.81) M 15.99(15.76, 16.04) rf 15.24(15.01, 15.29) rf 19.29(19.06,19.44) 
X-ray + FUV 90(83,96) 6.44(6.37, 6.52) c 14.44(14.38, 14.51) rf 15.86(15.80, 15.93) d 15.23(15.17, 15.30)'' 19.12(19.06,19.19) 
NOTE. — In our absorption line model, the transition oscillation coefficient and the damping factor are adopted from Vern er et al. I 
d!996h for O vn and O vm, and from lMortonl d2003h for O vi. " A is oxygen abundance in unit of solar. * The model predicated 
values. c The logarithmic value of To at the Galactic plane. d These uncertainty ranges can not be directly calculated from the joint 
analysis, and the error ranges listed here are scaled to have the same ranges as that of Mi. ' Assuming that all O vi absorptions are 
associated with O vn-bearing gas. See text for detail. 



The obtained A^h and T values are nearly identical to those 
from fitting the absorption lines alone. In addition, we con- 
strain L to be ~ 2(0.5,4) kpc. Fig. |2] presents the 68%, 90% 
and 99% confidence contours of L versus T. However, this 
isothermal modeling predicts an emission line intensity ra- 
tio of /ovii //ovm ~ 20, in contrast to the observed value of 
3 Q2.21 >. In the isothermal scenario, this ratio depends only 
on T (Eq|3l Fig. [3}. To reproduce this ratio, the tempera- 
ture needs to be ~ 10 6 33 K, which is significantly higher than 
what is inferred from the isothermal model (at > 99% con- 
fidence; Fig. 13. This inconsistency is apparently due to the 
over-simplification of the isothermality for the hot gas. 




6.10 6.15 6.20 6.25 

log[T(K)] 

FIG. 2. — The 68%, 90%, and 99% confidence contours of effective scale 
length (L) versus the temperature (T) of the emitting/absorbing gas at the CIE 
state, which is obtained by jointly-fitting the O VII (without O VIII) emission 
line intensity with the observed X-ray O VII and O VIII absorption lines. 

3.2. Non-isothermal Model 

We now consider a simple non-isothermal hot gaseous disk 
model to interpret the absorption and emission data. This 
model is motivated by the thick disk-like morphology of dif- 
fuse hot gas observed in nearby normal edge-on galaxies and 
inferred from our statistical sight-line dependence analysis of 
the Ne ix absorption line (Yao & Wang 2005; Tiillmann et 
al. 2006; see §4 for the discussion). In addition to the charac- 
terization of the density distribution as Eq. [2] the temperature 
distributions of the hot gas in the disk is assumed to be 

T = T e- z/(hTi) , (4) 
where the mid-plane value To and the scale height hj, along 
with h n and tiq in Eq. [2] are to be determined. This tem- 




5.5 6.0 6.5 7.0 7.5 8.0 

log[T(K)] 

FIG. 3. — The O VI, O VII, and O VIII line emissivities as a function of 
the gas temperature. The value for O VI has been scaled down by a factor of 
1000 for ease of demonstration. 

perature distribution is particularly motivated by the ap- 
parently lower temperatures of the absorbing gas toward 
LMC X-3 and Mrk 421 than those tow ard Galactic sources 
(lYao & Wangll2^05l:lwln7erai1l2^05h . 
From Eqs.|2]and|4] we can derive 

n = n (T/T r, (5) 

where j = hj /h n . Therefore, the differential hydrogen column 
density distribution is also a power law: 

cINh = ndL 

= ^L(T/T r~ l dT. (6) 

The corresponding ionic column density can be expressed as 

where A e and fi(T) are the element abundance and ionization 
fraction, and r m ; n is the minimum temperature considered in 
this work (see below). 

Similarly, the surface emission intensity is (Eq. Q]) 

A e f To , dEM 

where 

dEM _ L2ZV§7 (TV 1 - 1 
dT TqL \ToJ 

The intensity can then be calculated for an individual emission 
line or a broad band with the appropriately chosen A. 
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We modify our absline model to accommodate the power 
law temperature dependence of the absorbing hot gas column 
density (Eq. [7]l and also revise the emission line model con- 
structed in SjXT] according to Eq.[8] We then use these models 
to jointly fit the X-ray absorption and emission data (includ- 
ing both O vn and O vm emission lines). We consider the 
hot gas in the temperature range from T m \ n = 10 5 K to To; 
gas at lower (< 10 5 K) temperatures contributes little to the 
O vi, O vn, and O vm line absorption/emission. The absorp- 
tion/emission line joint fit with this non-isothermal model is 
as good as the one to the absorption lines alone Q3.lt . We ob- 
tain log[7b(K)] = 6.44(6.36, 6,54), L= 1.8(0.5,4.5) kpc, and 
7 < 1 .4. The fitted Nh and Vj values, as included in Table 12 
are consistent with those obtained in the isothermal modeling 
(§ 13.11 ). But the model predicated O vii/O vm emission line 
ratio, 3.7(2.2, 6.7), is now consistent with the observed value. 

It is constructive to specify how the individual parameters 
are constrained, although they are generally correlated in the 
spectral fitting. The relative O vn Ka to K/3 absorption line 
strength provides constraints mainly on Vb, Nova/Nowm on 7, 
Afovn on A^h (via the assumed Ao), the emission line intensity 
ratio /ovii/^oviii on 7b, and the ratio of Tovn to Abvn on L. 
Fig. Eta)-(c) presents the confidence contours of A^h, 7b, and 
L versus 7. It is clear that although the lower limit to 7 is 
poorly determined here, the L is well constrained to be < 10 
kpc at the 99% confident level. 

While the X-ray data do not give a tight constraint on 7, 
the inclusion of the FUSE O vi absorption observation in 
our fit can be helpful. In fact, the above non-isothermal 
modeling with the X-ray data alone over-predicts the O vi 
column density by a factor of ~ 2, judged by the best- 
fitted non-isothermal model (Table [2] j32.3t . If the O vi- 
bearing gas is entirely located in the hot gaseous disk, as 
prescribed for the X-ray absorption and emission, we can 
then obtain 7 = 0.6(0.4,0.9), L = 1.8(1.0,3.2) kpc, while the 
other parameters are included in Table [2] for comparison. 
The disk model parameters are log[7b(K)] = 6.44(6.37,6.52), 
n = 2.4(1.6,3.7) x 10" 3 cm" 3 , h n = 1.6(0.9, 2.9)/^ kpc, and 
hj = 1.0(0.6, 1.8)/£ kpc. The confidence contours for A^h, To, 
and L versus 7 are presented in Fig. |4] (d)-(f). Compared to 
Fig.[4](a)-(c), the tightening is partly due to the directly mea- 
sured \>b value for the O vi line. Of course, attributing the 
entire O vi absorption line to the hot gaseous disk represents 
an extreme case; the relationship between the O vi- and O vn- 
bearing gases is still very uncertain (see discussion in §4). 
But, if only part of the O vi absorption is associated with the 
disk, the lower limit to the 7 value will be increased further. 

As a consistency check, we further compare the predicted 
background emission intensity (Isxb) from the characterized 
exponential disk with the observed value ( 32. 2t . We first re- 
vise the thermal emitting gas model cemekl (as in XSPEC) 
and rename it as absmkl, 12 which takes To, Nh, L, and 
7 as input parameters and accommodates a power law EM 
of the gas, as presented in Eqs. |8]|9] Taking L = 2 kpc, 
7 = 0.6, and log[A r H (cm 2 )] = 19.12 (Table 0, we integrated 
this model with RASS response matrix to estimate the ex- 
pected hot gas emission intensity in the | keV band. Since 
the emission measure is mostly sensitive to the plasma tem- 
perature, we include the uncertainty in the mid-plane tem- 
perature, log[7b(K)] = 6.44(6.37,6.52) (Table©. We obtain 
/sxb = 1.0(0.4, 1.7) x 10" 4 counts s" 1 armin , which, within 
the errors, is consistent with the observed value ( 82. 2b . 

12 This model can be obtained from the authors on request. 



4. DISCUSSION 

4.1. Comparison with previous works 

We have presented a joint analysis of the z ~ absorption 
and emission lines in the X-ray spectra of Mrk 421 and found 
that the X-ray absorbing/emitting gas is consistent with being 
Galactic, most likely located in a thick hot gaseous disk with 
an effective scale length ~ 2 kpc. This same disk may also 
explain much of the O vi-bearing gas observed in the far-UV 

vi line absorption. Here we compare our results and con- 
clusions with those from existing studies on the hot gas along 
the same sight-line; some of which also give estimates of the 
scale length of the X-ray absorbing/emitting gas. All existing 
studies assume that the hot gas along the sight-line is in an 
isothermal state. We have shown that this assumption is prob- 
lematic, at least in the Galactic interpretation. Instead, our 
analysis have allowed for nonthermality, which is facilitated 
by the joint-fit approach. This approach maximizes the use of 
the information in the data, and propagates the error of each 
model parameter and takes care of the correlations between 
different parameters automatically, and thus yields more reli- 
able results. There are, of course, other differences in both the 
analysis and interpretation, sometimes leading to very differ- 
e nt conclusions in v a rious e xisting studies. 

iRasmussen et al. I (120031) . for example, estimate a scale 
length > 140 kpc, based on the averaged O vn and O vm 
line absorptions along the sight-lines toward three AGNs 
(Mrk 421 3C 273, and PKS 2 155-304) and the O vn emission 
line from McCamm on et al. 1 d2002f) . This estimate uses an 
ionization fraction fi of 0.5 for O vn and an oxygen abundance 
of 0.3 solar. The large scale is in contrast to our inferred value 
L ~ 2 kpc (§3). The bulk of this discrepancy can be easily ac- 
counted for (see Eq. 3): oxygen abundance (0.3 solar versus 

1 solar), hydrogen density distribution (uniform versus expo- 
nential, causing a factor of 2 difference), and the mean hot 
gas temperature (6.17 versus 6.35 in logarithm, giving a fac- 
tor of ~ 6 due to the change of both ionic fraction and emis- 
sivity). These differences give a combined factor of ~ 40; 
the different adopted atom data may also contribute. The 
higher average gas temperature adopted in Rasm ussen et al. I 
(2003) results from the inclusion of the 3C 273 sight-line that 
passes through the inner region of the Galaxy and is appar- 
e ntly contaminated by th e North Polar Spur or the radio Loop 
I dSnowden et al. l[T997l) . Additional difference likely arises 
fr om the non-isothermality of our hot gaseous disk model 
( 33.2t . An absorption line tends to be produced by gas at 
lower temperatures than the corresponding emission line of 
the same transition (e.g., O vn Ka; compare Fig.[3]and Fig. 1 
in Yao & Wang 2005). 

Williams et al. (2005) have also compared the z ~ X-ray 
and far-UV (O vi) absorption lines toward Mrk 421, under the 
isothermal gas assumption. They find that the velocity disper- 
sion inferred from the X-ray absorption lines is smaller than 
that of the direct measurement of the O vi line. As an alterna- 
tive, an extragalactic origin of the X-ray-absorbing gas with 
an inclusion of the photo-ionization process is explored and 
is found to be consistent with the absorption data, which re- 
quires a super-solar abundance ratio of Ne/O. In this work, we 
only consider the oxygen absorption lines to avoid the uncer- 
tainty in the ratio. We find that, while our measured equivalent 
widths (EW s) for other oxygen line s are consistent with those 
reported by IWilliams et al. I (120051) . the EW of O vn Ka line 
[11.4(10.20, 12.68) mA] is about ~ 20% hig her. Note that 
our value is consistent with that obtained by Kaastr aet al. I 
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Fig. 4. — The 68%, 90%, and 99% confidence contours of the non-isothermal model parameters. They are obtained from the joint analysis 
the X-ray absorption/emission data with (right panels) and without (left panels) the inclusion of the far-UV O VI absorption line. See text for 
details. 
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(2006). We believe that this discrepancy is mainly due to their 
improper weighted instrumental response file 13 for the final 
co-added spectrum; the higher (> 1) order confusion of the 
HRC observation (if not properly being taken into account), 
the different placement of the continuum, and the different 
software packages utilized in the data analysis may also con- 
tribute. A higher O vn Ka/K/3 EW ratio indicates a small 
saturation of the O vn Ka absorption and therefore requires 
a large b v , which then becomes consistent with the value of 
the O vi line. Jointly analyzing the absorption and emission 
data, we have shown here that the data can be explained by 
the Galactic hot gas. 

Savage et al. (2005) have explored the FUSE spectra 
of Mrk 421 and two nearby stars BD +38°2182 (l,b = 
182.°16,62.°21) and HD 93521 {l,b = 183.° 14, 62.° 15) in great 
detail. They have examined the relationships of the absorption 
by O vi to those by lowly-ionized species (e.g., H i, C n, C in, 
and N n) as well as to those by highly-ionized ones (O vn and 
O vm). The observed "broad" O vi absorption is divided into 
the two different velocity components, from -140 to 60 km 
s" 1 and from 60 to 165 km s" 1 . A good correlation between 
the positive high-velocity O vi component and low ionization 
gas absorptions is found. Interestingly, the high-velocity O vi 
component does not seem to show up in the spectra of the two 
stars. If these sight-lines are not atypical, the positive high 
velocity absorption should then occur at a distance at least 3.5 
kpc off the Galactic plane. In this case, the absorption cannot 
arise from the same location as the O vn- and O vm-absorbing 
gas. 

Some portion of the observed O vi absorptions should 
be associated with the O vn-bearing gas, which we as- 
sume to be located in a thick hot gaseous disk. To fur- 
ther probe this association, we reprocessed the FUSE ob- 
servations of BD +38°2182, following the standard CalFuse 
scripts. We find that the main O vi absorption along the sigh 
line can be characterized as a narrow component with Vb = 
33(3 1 , 35) km s" 1 and log[jV OV i(cirr 2 )] = 14.20(14. 17 , 14.22). 
However, because the continuum placement is very arbitrary 
(see Zsargo et al. 2003 and Savage et al. 2005 for further dis- 
cussions), this characterization is very uncertain. A careful 
check suggests that up to 30% of the broad O vi line absorp- 
tion observed in Mrk 421 can be superposed on the narrow 
line easily seen in the spectrum of BD +38°2182. This "com- 
mon" broad component could be associated with the O vn- 
bearing gas. The remaining O vi line absorptions in the sight- 
lines toward the star and of Mrk 42 1 could then arise in the 
discrete interfaces at bounda ries of cool clouds, emersed in 
hot gas ( Savag eet al. 1 l200"5l) ; the apparent broadening ob- 
served in Mrk 42 1 may be due to the superposition of different 
velocity components over a large distance. The true absorp- 
tion and emission of the gas should be a convolution of the lo- 
cal temperature and density dispersions with the global depen- 
dences on the off-plane distance. But currently there is little 
observational or theoretical guidances as to the specific forms 
of the dispersions. Our simple exponential model, hopefully, 
provides a first-order correction to the isothermal assumption, 
or at least demonstrates that the non-isothermality of the hot 
gas should be seriously explored. 

To examine the influence of the "common" component on 
our results, we perform the same joint analysis as in 33.21 but 

13 The response files should be weighted with the corresponding count 
spectrum rather than with the exposure time only, since the flux of the source 
is varying significantly. 



requiring the hot gaseous disk accounts for only 30% of the 
total O vi absorption detected in the Mkr 421 sight-line. The 
constrained parameter ranges are very similar to those illus- 
trated in Fig. |4] (d)-(f), except for a slightly larger 7 and a 
smaller To values. Clearly, we cannot unambiguously deter- 
mine the exact fraction of the observed O vi absorption that 
may be associated with X-ray absorbing gas. But this uncer- 
tainty is not expected to qualitatively affect our results and 
conclusions about the overall non-isothermality and extent of 
the X-ray-emitting/absorbing gas. 

4.2. Location of the O vn-bearing gas 

The above analysis and results show that the gas that pro- 
duces the X-ray line absorption and emission along the sight- 
line toward Mrk 421 is consistent with a Galactic hot thick 
disk interpretation. But based on the data for this single sight- 
line, we cannot conclusively determine the location of the gas. 
As mentioned above, alternative scenarios (i.e. Galactic halo 
and Local Group) of the z~0O vn-bearing gas have been 
proposed, allowed for by the limited spectral resolution of the 
X-ray grating observations. So let us review and discuss the 
merit and implic ations of these p r oposa ls. 

For instance, iNicastro et al. I (|2002) and Williams et al. 
(120051) adopted a Local Group WHIM interpretation and 
found a consistent scale length of ~ 1 Mpc of the gas in the 
sight lines of PKS 2155-304 and Mrk 421. This consistency 
is partly a natural outcome of the low-density assumption, 
hence the importance of the photo-ionization by the extra- 
galactic background in determining the ionic fractions of the 
gas. This interpretation, however, requires the WHIM to be 
either metal-rich (e.g., Ao > 0.3 solar) or severely (spherical) 
asymmetric (i.e . , a small sky-covering- fraction f c ), or both 
(Williams et al. 2005; Fang et al. 2006), because otherwise 
the derived total baryonic matter in the Local Group would 
be even higher than its gravitational mass (~ 2 x 10 12 Mp,; 
Courte au & van den Bergh 1999). Recently, Fan g et aL~l 
(2006) performed a systematic search for the O vn absorp- 
tion line in AGN spectra from archived Chandra and XMM 
data and found that the O vn absorption line is detected sig- 
nificantly as long as the counting statistic is reasonably good. 
They concluded that f c is >63% and likely to be unit. This 
large f c is clearly contradictory to the requirement of the 
WHIM interpretation that f c must be <30% even a ssuming 
the solar metallicity for the hot gas (Fang et al. 2006). There- 
fore, the WHIM scenario is very unlikely, if not impossible, 
to explain the bulk of the O vn-bearing gas. 

Alternatively, hot gas is also expected to be present in the 
Galactic halo on a scale of ~ 100 kpc. Such a hot gaseous 
Galactic halo is predicted in galaxy formation theories. Such 
gas is heated at the so-called virial shocks and through sub- 
sequent gravitational compression (e.g., Birnboim & Dekel 
2003). This gas, cooling radiatively, can maintain lasting star 
formation in galactic disks. But simulations tend to over- 
predict the cooled gas content in galactic disks, which is the 
so-called "over-cooling" problem (e.g., Navarro & Steinmetz 
1997). The inclusion of galactic energy feedback may alle- 
viate the problem, although this has not been satisfactorily 
demonstrated. Observationally, the existence of a hot gaseous 
halo around our Galaxy is supported by observations of high 
velocity (HV; |v L sr| > 100 km s" 1 ) clouds (HVCs). Some of 
these clouds, especially compact ones, appear to be located 
far away from the Galactic plane, although tight distance con- 
straints are very few. The hot gaseous halo may be needed to 
confine such clouds and may even be their origin through con- 
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FIG. 5. — The expected surface intensity distribution in 0.3-1.5 keV band 
as a function of the vertical distance off the Galactic plane. 

densation. Furthermore, HV O vi absorption lines are some- 
times associated with HVCs and are bel ieved to be produce d 
at their interfaces with the hot halo gas (Sem bach et al. 120031) . 
So far, however, there is little observational evidence for a 
large-scale diffuse X-ray-emitting or -absorbing halo around 
our Galaxy or other similar disk galaxies (see Wang 2006 for a 
review), consistent with the conclusion reached in the present 
work. A scenario that reconciles these observations is that 
much of the hot halo gas is very low in metallicity, because 
both the X-ray emission and absorption arise from metal ions. 
Such a Galactic hot gaseous halo with low metallicity may be 
a natural product of the strong selective radiative cooling if the 
accreted gas has a distribution of metal abundances; the gas 
with a higher metal abundance cools faster and forms denser 
clouds (probably HVCs), whereas the remaining gas tends to 
have zero or low metallicity, especially in inner regions of 
galactic halos. This scenario explains the low X-ray emis- 
sion/absorption of the hot halo gas and further alleviates the 
over-cooling problem. A detailed modeling of the scenario 
will be presented elsewhere. 

We favor the interpretation that the bulk of the z~0O vn- 
bearing gas represents a hot thick gaseous disk around the 
Milky Way, consistent with previous studies as well as the 
present one. In this case, the gas represents the hot phase 
of the ISM and is presumably heated primarily by supernova 
blastwaves. The presence of the gas is also traced by O vi 
absorption. As its ionic fraction peaks at 3 x 10 5 K, O vi 
is believed to present primarily at interfaces between cool 
and hot gases and/or in hot gas cooling. The O vi line ab- 
sorption is observed ubiquitously at relatively low velocity 
(LV; |vlsr| $100 km s" 1 ), clearly indicati ng a predominantly 
Galactic disk origin (IWakker et al. I I2003T) . The spatial dis- 
tribution of the O vi absorption indeed has a characteristic 
vertical scale height of ~ 2.3 kpc around the Galactic plane 
( Sava ge et al. 120 03) and is thus in a good agreement with that 



of the hotter O vn-bearing gas. Therefore, the presence of the 
hot thick gaseous disk is fully expected and provides a natu- 
ral explanation for much of the X-ray/far-UV absorption and 
emission. 

The hot thick gaseous disk interpretation is also supported 
by the X-ray emission observations of nearby disk galaxies 
like our own. Both the intensity and spatial extent of diffuse 
X-ray emission appear to be correlated with star formation in 
galactic disks. Typically the emission is confined to galactic 
coronae arou nd the disks with vertic al extent no more than a 
few kpc (e.g jTuilmann et al. I (120061) . Wang et al. 2003). For 
ease of comparison with these studies, we present in Fig. [5] 
the emission intensity distribution as a function of the ver- 
tical distance away from the Galactic plane with the best-fit 
parameters of the hot thick gaseous disk ( {33.21 i. Because the 
7 is not well constrained, the distribution for several specific 
7 values are included. For a very broad range of 7, the emis- 
sion is essentially concentrated within a few kpc off -plane dis- 
tance. This concentration is consistent with those findings in 
the ne arby disk galaxies (Wa ng et al. Il2003t iTullmann et al. I 
2006), suggesting a common behavior and a common nature 
of hot gas in these galaxies. 

To summarize, we have carried out a systematic inves- 
tigation of the available X-ray O vn and O vm absorp- 
tion/emission tracers of Galactic diffuse hot gas in the field 
of Mrk 421. A joint analysis of the data, complemented by 
the FUSE O vi measurement, enables us for the first time 
to self-consistently and simultaneously characterize the non- 
isothermality and vertical scale height of the hot gas in the 
sight-line toward Mrk 421. We have established that the X- 
ray-emitting/absorbing hot gas cannot be in an isothermal 
state. But the X-ray data can be explained with a hot gaseous 
disk model, which is meant to provide a simple, yet physi- 
cally plausible characterization of the global temperature and 
density distributions of the hot gas. Indeed, our inferred scale 
height of the hot gaseous disk and the gas density at the Galac- 
tic plane are consistent with our previous estimations based 
on the Galactic latitu de dependence of G alactic Ne ix line ab- 
sorptions observed dYao & Wang | |2005). 
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